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Abstract 

The 2014-2015 Ebola Virus Disease (EVD) epidemic in West Africa was the largest ever recorded, 
representing a fundamental shift in Ebola epidemiology with unprecedented spatiotemporal complex¬ 
ity. We developed spatial transmission models using a gravity-model framework to explain spatiotem¬ 
poral dynamics of EVD in West Africa at both the national and district-level scales, and to compare 
effectiveness of local interventions (e.g. local quarantine) and long-range interventions (e.g. border- 
closures). Incorporating spatial interactions, the gravity model successfully captures the multiple 
waves of epidemic growth observed in Guinea. Model simulations indicate that local-transmission re¬ 
ductions were most effective in Liberia, while long-range transmission was dominant in Sierra Leone. 
The model indicates the presence of spatial herd protection, wherein intervention in one region has 
a protective effect on surrounding regions. The district-level intervention analysis indicates the pres¬ 
ence of intervention-amplifying regions, which provide above-expected levels of reduction in cases 
and deaths beyond their borders. The gravity-modeling approach accurately captured the spatial 
spread patterns of EVD at both country and district levels, and helps to identify the most effective 
locales for intervention. This model structure and intervention analysis provides information that 
can be used by public health policymakers to assist planning and response efforts for future epidemics. 
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1 Introduction 

The outbreak of Ebola virus disease (EVD) in West Africa caused 28,646 cases and 11,323 deaths 
as of March 30 , 2016 [ ]. The current outbreak is of the EBOV (Zaire Ebola Virus) strain, the most 
fatal strain [2, 3]. The largest previous outbreaks of the EBOV strain occurred in the Democratic 
Republic of the Congo in 1976 and 1995 , causing 318 and 315 cases [4]. 

Clinical progression includes two broad stages of infection, often characterized as early and late 
[5]. In the first stage, approximately five to seven days, symptoms include fever, weakness, headache, 
muscle/joint pain, diarrhea, and nausea [5, 6]. In some patients, the disease progresses to a second 
stage, with symptoms including hemorrhaging, neurological symptoms, tachypnea, hiccups, and 
anuria [5, 7]. Mortality rates are higher among those exhibiting second-stage symptoms [5, 7]. 
EVD is transmitted through direct contact with an infected individual [8]. Transmission risk 
factors include contact with bodily fluids, close contact with a patient, needle reuse, and contact 
with cadavers, often prepared for burial by the family [8—10]. 

This outbreak was primarily in the contiguous countries of Guinea, Sierra Leone, and Liberia, 
which experienced widespread and intense transmission [1, 11, 12], Interventions have included 
quarantine, case isolation, additional treatment centers, border closures, and lockdowns, restricting 
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Figure 1: A diagram of the model structure (left panel) and a map of West Africa showing the locations 
of each capital used as a population center for each patch (right panel). In each country, the population is 
compartmentalized into the following categories: susceptible (5), exposed (E), infected in stage one (Ji), 
infected in stage two (I 2 ), died but not yet buried (F), and recovered ( R ). Transmission is possible from 
persons in compartments I±, I 2 , or F of any country, to persons in compartment 5, either within the country 
or to the other countries. 


travel within a region (as in a military-enforced curfew) [13]. Public health authorities must allocate 
resources effectively, focusing personnel and funds to respond best to outbreaks. However, not all 
response measures are equally beneficial or cost-effective, and testing the relative benefits of each is 
often impossible or unethical. Modeling thus provides a valuable tool for comparing interventions 
and identifying areas where interventions are most effective. Modeling analyses can inform public 
health policy regarding ongoing and future outbreaks. 

Dynamic spatial modeling has been proposed as a useful approach to understand the spread of 
EVD and evaluate response strategies’ benefits [ ]. The recent outbreak has sparked an increase 

in EVD modeling [15-18]. The spread between contiguous countries in the current EVD outbreak 
highlights the spatial element to its proliferation [ ]; however, as previous EVD outbreaks were 
more localized than the 2014-2015 epidemic [19], there is little historical data on the geospatial 
spread of Ebola. Mobility data, which may help inform spatial EVD modeling, is limited, although 
some studies have highlighted its usefulness and extrapolated based on mobility data from other 
regions [14, 20]. 

To address the need for spatial EVD models, models have been developed to examine local spa¬ 
tial spread within Liberia [ 14,1 , ] and evaluate the potential risk of international spread using 

data such as airline traffic patterns [21, 23, 24]. Our models add to existing spatial models by in¬ 
corporating intervention comparison and an exploration of the dynamics between different regions, 
using an uncomplicated model structure. Between-country mobility is important in this epidemic 
because the borders between countries are porous; borders are drawn across community or tribal 
lines, resulting in frequent border-crossing to visit family, conduct trade, or settle disputes [25]. 
In addition, Santermans and coauthors recently demonstrated that the outbreak is heterogeneous: 
transmission rates differ between locales [2 ]. Our district-level model offers a mechanistic expla¬ 
nation of observed spatial heterogeneities, adding to the existing literature by explicitly simulating 
the interactions between districts. 
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In this study, we present spatial models of EVD transmission in West Africa, using a gravity 
model framework, which captures dynamics of local (within-region) and long-range (inter-region) 
transmission. Gravity models are used in many population-mobility applications, especially to 
relate spatial spread of a disease to regional population sizes and distances between population 
centers [27]. Gravity models have been applied to other diseases, including influenza in the US 
and cholera in Haiti [27, 28], and have been used to examine general mobility patterns in West 
Africa[l I]. 

We demonstrate that the gravity modeling approach fits and forecasts case and death trends in 
each country, and describes transmission between countries. The district-level model successfully 
simulates local geospatial spread of EVD, indicating that the gravity model framework captures 
much of the local spatial heterogeneity. Using the model structures, we examined interventions at 
country and district scales, evaluating the relative success of intervention types as well as the most 
responsive regions. 


2 Methods 

2.1 Country-level Model Structure 

We developed a compartmental gravity model using ordinary differential equations to model spa- 
tiotemporal progression of EVD (Figure 1, model equations and details given in the Supplemen¬ 
tary Information). Each country’s compartmental structure includes susceptible (S), latent ( E ), 
two-stage infection (ii, I 2 ), funeral (F), recovered (R), and deceased (D), based on previous com¬ 
partmental models [29, 30]. Patients in stage I\ can recover or transition to 1-2, where they may 
recover with lower probability or transition to F. based on clinical observations of symptom pro¬ 
gression and mortality [7]. The latent period and two-stage infection are based on the clinical 
progression of EVD, wherein patients become contagious in I \. with increasing contagiousness in I 2 
[5, 8, 15, 30, 3 ]. Funerals play a role in transmission because of the high viral load of the deceased 
and frequent contact with susceptible persons due to cultural burial practices [32]. The precise 
relative magnitude of this contribution is unknown; previous models have shown that the relative 
contributions of each stage are unidentifiable from early data [30, 33]. Thus, patients in F were 
assumed to be as contagious as those in F- 

The spatial component of the model is a three-patch gravity model. Each country is one patch, 
with a compartmental model within it and the capital as the population center, similar to previous 
gravity models, which used political capitals/large cities as centers [27, 28]. After EVD appeared 
in western Guinea, subsequent cases were in the capital, Conakry, suggesting that the capital acts 
as the central population hub. Similar progression occurred in each country: after EVD cases 
appeared in border regions, cases soon appeared in the capital [34], 

The force-of-infection term for patch i consists of transmission from within the region and 
transmission from other patches into patch i. Each long-range transmission term is determined 
by a gravity term, proportional to the population sizes, and inversely proportional to the squared 
distance between them. 

Cumulative local cases were measured by the number of infections in a patch due to local 
transmission. Cumulative long-range cases were measured by number of cases due to long-range 
transmission. 
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2.2 District-level spatial model 

We developed a more detailed gravity model of West Africa incorporating the 14 districts/areas of 
Sierra Leone, 15 counties of Liberia, and 34 prefectures of Guinea. This model is structured the 
same way as the country-level model: a compartmental model in each administrative unit, linked 
to all other patches by gravity terms. Transmission is separated into local (within-district) and 
long-range (between-district). 

2.3 Data and Simulation Setup 

Case and death incidence in each country was collected from World Health Organization (WHO) 
situation reports on EVD from March 29 to October 31 , 2014 [ ]. The country-model simulations 
used data from May 24 through September 30 for fitting model parameters. Data from October 
was used for validation, to compare model projections to data unused in fitting. Road distance 
between centers was used for the gravity component [35]. We evaluated direct distance between 
capitals in the country-level model, yielding similar results. 

For further verification of the model, we also tested the model’s ability to fit the outbreak when 
data was incorporated from March 29. The fits and forecasts from these simulations were similar 
to those incorporating data from May onward (Supplemental Figures 1-2). 

The district model was simulated from March 31 , 2014 to January 31 , 2015 . The model was 
compared to data on the outbreaks local geospatial progression from each district, using incidence 
of cases and deaths from WHO updates and the UN OCHA database [34, 36]. 

2.4 Parameter Values and Estimation Methods 

Model parameters were determined from clinical literature and fitting to available data, similarly 
to previous models [18, 30, 37]. In the country-level model, nine parameters were fitted to the data: 
transmission (/3) and death (A) as well as three gravity-term constants, n, were separately fitted 
for each country. . The parameter (/c) adjusts the gravity terms to reflect the balance of local and 
long-range transmission in each country. Parameters, definitions, units, ranges, and sources are 
given in Supplementary Table SI. 

In the country-level model, 500 sets of initial values for all parameters were selected with Latin 
Hypercube (LH) sampling from realistic ranges for parameters, determined from WHO, CDC, and 
literature data, given in Supp. Table SI [5, 7, 16, 29, 30]. For each of the 500 parameter sets, only 
transmission rates, death rates, and k were fitted by least-squares, using Nelder-Mead optimization 
in MATLAB [38]. All other parameters were held constant to the values from the LH sample. 

In the district-level model, best-fit parameters from the country-level model were used for all 
parameters except k norm . A different k norm was fitted for each district to reflect differing at-risk 
populations and reporting rates. Since fitting 63 parameters is highly complex, slow, and possibly 
unidentifiable, k norrn parameters were calculated according to final size, then adjusted by hand to 
generate curves matching the progression of the outbreak, both overall and in each district. While 
all other initial conditions were determined from early data and k norm , in Conakry, Guinea, the 
initial conditions were fitted independently from k n0 rm , reflecting possible reporting rate differences 
between initial situation reports and subsequent data. 

2.5 Intervention Simulations 

Interventions were compared by the number of cumulative cases the model forecasted on October 
31 for different types and levels of intervention. Reduction of local transmission represents inter- 
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Figure 2: Cumulative cases and deaths for each country. Data depicted with black circles (from May 24, 
2014 to September 30, 2014) was used for model fitting to the data. Data depicted with red crosses (from 
October 1 to October 29) is displayed for validation of the model fits. The model projections are from May 
24 to October 31, 2014. The black circles represent the data used for fitting; the red crosses represent the 
data used for validation. The line represents the models overall best fit; the shaded regions represent the 
95% confidence interval. 


ventions that limit contact of an EVD patient with others in his or her home country, including 
quarantine, isolation, hospitalization, and local lockdowns. Reduction of long-range transmission 
represents border closures, large-scale lockdowns, or other intervention measures that reduce the 
transmission between countries. 

In the district-level model, interventions were simulated by eliminating the outbreak in one 
patch (both local and long-range intervention), and simulating to January 31. This is analogous to 
effective case isolation and quarantine measures within a district, which would limit the transmission 
of EVD within and from that district. The effectiveness of intervention in one district was measured 
by calculating the cumulative number of cases reduced in all other districts. The percent reduction 
from intervention in a district and percent reduction relative to size of the outbreak in that district 
were calculated. 


3 Results 

3.1 Country-level Model Fits and Predictions 

The best-fit accurately fit and forecasted the outbreak data and trends for each country (Figure 
2). The best-fit model forecast of cases on October 29, 2014, was 14,070; the WHO case data on 
that date was 13,540 [39]. 
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Figure 3: Cases from local transmission and long-range transmission. Local transmission, from transmission 
within the country, is depicted in the color of the country in Figure 1. Long-range transmission, or cases 
from outside the country, is depicted in gray. Model projections are from May 24 to October 31. 


3.2 Local versus Long-Range Transmission 

Cumulative cases in each country due to local and long-range transmission are shown in Figure 3. 
Liberia had more local than long-range transmission. Early on, Guinea had more local transmission, 
due to initial local cases. As the epidemic progressed, the transmission ranges overlapped, although 
the best-fit trajectory for long-range transmission remained smaller than the local-transmission 
contribution. In Sierra Leone, the best-fit trajectory of long-range transmission was more significant 
than local transmission, although forecasted ranges were similar. 

3.3 District-level Model 

The district model was successful in fitting outbreak data for cases and deaths in each district. The 
model captured the final size of the outbreak in each patch with an R 2 value of 0.96 (Figure 4) 
as well as matching the progression of the outbreak (Figure 5). The model captured the temporal 
dynamics of districts with negligible case counts or high levels of transmission. The time-course 
plots of all 63 district-level fits are included in the Supplement, as is a scatter plot of all data vs 
model values ( R 2 = 0.83). 

The district-level model accurately captured the trend of spread of EVD through second-level 
administrative units. While overall patterns at the district level are captured and show the correct 
ordering, the actual speed of disease spread in the model was faster in some districts than in the 
data, likely due to a combination of stochastic introductions and reporting delays. The district-level 
model was able to forecast local spread of EVD, matching data on spatial progression of EVD to 
different locales as the outbreak intensified (Figure 5). An animation of the district-level simulation 
is given in the Web Supplement (Video SI). 

3.4 Intervention Simulations: Country-Level 

According to our intervention simulations (Figure 6), reduction in local transmission, from inter¬ 
ventions such as isolation or improved case-finding and Ebola treatment unit (ETU) capacity, is 
most effective in Liberia. In the model, eliminating local transmission in Liberia reduced the out¬ 
break by up to 11,000 cases in all countries by October 31, 2014, a 76% reduction. Reduction of 
long-range transmission was most effective in Sierra Leone. Eliminating long-range transmission 
into Sierra Leone reduced the outbreak by up to 9,600 cases in all countries by October 31, a 66% 
reduction. 
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Figure 4: The final size of the outbreak in each of the 63 patches: comparison of model values (x axis) and 
data (y axis). R-squared value of 0.9635. 

3.5 Intervention Simulations: District-Level 

Intervention analysis identified regions in which intervention was most successful (Figure 7). Over¬ 
all, the most successful regions for intervention were Conakry, Coyah, and Dubreka, Guinea, and 
Lofa, and Monrovia/Montserrado, Liberia. 

To evaluate whether intervention effects were disproportionately significant compared to cases 
within the target region, we defined the Intervention-Amplifying Region (IAR) level by the ratio 
of cases prevented in other patches to cases within the patch, for regions with more than 0.05% of 
the epidemics cases. Dubreka, Coyah, and Conakry were the most effective IARs (Figure 7). 

4 Discussion 

Our results demonstrate that outbreak dynamics in all three countries can be accurately captured 
using a gravity-model approach. Moreover, country-level model forecasts successfully predicted 
cumulative cases and deaths one month ahead (Figure 3). The models were able to capture simul¬ 
taneously the dynamic interactions within and between each regionat both country and district 
levels. 

4.1 Spatial herd protection: intervention in one region protects others 

The model simulations suggest a form of spatial herd protection, wherein interventions in one region 
benefit surrounding regions as well, by reducing the spatial transmission between them. 

In the country-level model, we compared the effects of local and long-range interventions in each 
country. Local interventions include more strict quarantine and isolation procedures, increased case¬ 
finding and ETU capacity, safer burial practices, and other interventions that reduce contact of 
susceptible persons with infected persons in the local community, as well as behavioral changes 
that reduce local transmission. The model indicates that these interventions were most effective 
in Liberia, both in reducing the outbreak in Liberia and mitigating the whole epidemic. Indeed, 
Liberias outbreak had the fastest initial growth rate, it was the first to turn over and end [ ]. Our 
results suggest that Liberian interventions may have had significant indirect effects on the epidemic 
dynamics in Guinea and Sierra Leone as well. 


7 




Data 



Model Data 



Model 





I- 11 Cases 

II- 100 Cases 
101-200 Cases 
201-300 Cases 
301+ Cases 
No Data 



Figure 5: Map of geospatial spread of Ebola as determined by the district model (right) and WHO data 
(left). The shaded regions are the regions with cases of EVD; the intensity represents intensity of transmission 
(number of cases). Gray regions represent districts for which no data is reported on the given date. On the 
May 26/28 figure, May 26 is used for Sierra Leone and May 28 is used for Guinea. 



Case Reductions for Different Levels of Transmission Reduction 



Figure 6: Cases reduced by various interventions at the country-level. Transmission reduction represents 
reduction of local or long-range transmission in increments of one percent. Case reduction is measured 
by comparing total cases in all 3 countries in the presence of intervention to total cases in the absence of 
intervention, on October 31, 2014. 


The most effective long-range transmission reduction was in Sierra Leone, with 66% case reduc¬ 
tion across all countries when long-range transmission was completely eliminated. This reflects the 
porous border in Sierra Leone. Based on the district-level intervention analysis, no single district 
in Sierra Leone acts as a source of long-range transmission to other locales, rather all districts 
play some combined role. Thus, the country-level impact of long-range transmission is possibly 
a summative combination of two factors across all districts: long-range transmission into Sierra 
Leone and early introduction of cases from Sierra Leone into Liberia. 

In the district-level model, the capitals of Montserrado, Liberia and Conakry, Guinea were 
highly effective intervention sites. Their large population size makes them important in dictating 
the dynamics of the outbreak in surrounding regions. Dubreka and Coyah, Guinea, near Conakry, 
were effective intervention sites despite low numbers of cases. This indicates that they influenced 
the larger outbreaks surrounding them, especially Conakry. Intervention in these districts was 
especially significant in reducing the scope of the outbreak, not just within the district but also in 
surrounding regions. 

4.2 Intervention-amplifying regions: outsized impact on the outbreak dynamics 

In some regions, intervention produced unexpectedly significant reductions in the overall outbreak: 
these districts amplify the effects of intervention to reduce disproportionately many cases outside 
the regionthese are denoted Intervention-Amplifying Regions (IARs). The most significant IARs 
include Conakry, Coyah, and Dubreka in Guinea and Montserrado in Liberia (Figure 7). For in¬ 
stance, Dubreka made up only 0 . 13 % of the outbreak but intervention in Dubreka alone reduced 
the size of the outbreak in all regions by 21%, drastically higher than its case contribution would 
suggest. Coyah and Dubreka were the most dramatic IARs because they had relatively small num¬ 
bers of cases. Their proximity to Conakry, which had a large number of cases and also functioned 
as an IAR, implies that Coyah and Dubreka could have influenced outbreak dynamics in Conakry. 

The outsize importance of these districts in the dynamics of the current outbreak makes them 
a possible target for increased intervention during future outbreaks in West Africa. 

4.3 Model fits unique case curve in Guinea 

Several papers have noted that Guineas peculiar outbreak curve is difficult to fit using simple 
models, due to plateaus in cumulative incidence between growth periods [18, 24, 40]. However, 
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Figure 7: Cases reduced by various interventions at the district-level. Transmission reduction consists of 
elimination of outbreak in one district. Case reduction is measured by (a) percent reduction of outbreak 
in other 62 patches and (b) percent reduction in other 62 patches relative to size of the outbreak in the 
intervention site (Intervention-Amplifying Region score). Panel (c) displays the IAR scores on a scatter plot, 
with the associated quartiles depicted by the shaded regions. The blank region represents the patches with 
insufficient cases to be considered for intervention. 


by including spatial interaction between countries, we successfully captured Guineas outbreak dy¬ 
namics (Figure 2). This suggests that spatial interactions, such as local die-outs of EVD followed 
by long-range reintroductions, may be responsible for the unusual incidence patterns in Guinea. 
The country-level models ability to capture the singular epidemic progression in Guinea, due to its 
spatial transmission component, demonstrates that long-range dynamics do not just affect intro¬ 
duction, but also later stages of the outbreak. 

4.4 Discussion of district-level model 

The district-level model, which separates the affected countries into 63 patches based on secondary 
administrative units, was successful in capturing the overall dynamics of spatial transmission. It 
captured the spread of EVD from the initial sites of the outbreak to other locales in the region 
in an order and magnitude similar to the actual progression of the epidemic. This suggests that 
the county-level model captures the underlying transmission patterns that led to spread of EVD 
between different locales in the affected countries (Figure 5). These results are perhaps surprising 
given that the district-level model uses the same transmission parameters for all patches within the 
same countrythis suggests that even though there are likely to be significant local heterogeneities 
from patch to patch [41], most of these variations can be captured using the relatively simple 
framework afforded by the gravity model. 

In fact, a model that incorporates local heterogeneities often requires a large amount of case 
data at the local scale [ ]. This gravity model provided a reasonable approximation of the early 

dynamics of the outbreak, based simply on initial conditions from March 30, soon after surveillance 
began: it can capture local differences in transmission without complex (difficult-to-simulate) mod¬ 
els or extensive datasets. Thus, it could be applicable for early-outbreak forecasting by providing 
an indication of the areas at greatest risk for Ebola cases, in order to guide resources and aid to 
those locations. 
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4.5 Simplifying Assumptions and Limitations 

Our models focus on early proliferation throughout West Africa: interactions between regions, 
which led to the spatial spread of Ebola. Regional EVD introductions are inherently stochastic; 
thus the model did not capture the precise temporal pattern in West Africa in every district. 
Stochastic model simulations might better characterize the randomness of regional introduction, 
but a stochastic implementation of the 63-patch regional structure would be computationally de¬ 
manding. Therefore, a deterministic model offers an elegant method to provide insight into the 
role of spatial dynamics. 

In addition, we considered the time scales that allowed us to capture the spatial dynamics of 
interest. In the country-level model, the turnover and re-ignition of the epidemic in Guinea was of 
particular interest, since simple compartmental models provide little insight into its mechanism [ ]. 

Our time-course included the first ignition, turnover, and re-ignition; since simpler compartmental 
models fail to capture the dynamics of the outbreak, our models highly accurate fit implies that 
the spatial structure captures the dynamics affecting Guinea. 

Our time-course ends in November-December 2014, at which time interventions were increasing 
throughout West Africa. This massive intervention scale-up affected the parameters for the out¬ 
break: in order to simulate the progression of the outbreak after this time, time-variant parameters 
could be used. The late dynamics were not the focus of this research; however, the model could be 
implemented with time-variant parameters, potentially providing insight into late dynamics of the 
outbreak. 

4.6 Future research directions 

While this model captures epidemic dynamics without incorporating explicit movement patterns, 
further mobility data (e.g. from cell phone carriers) could elucidate dominant movement patterns 
in West Africa. This data could provide a validation of the gravity models success, if the gravity 
model accounts for these patterns. 

The model does not account for a detailed population structure: an organization of individ¬ 
uals into households, villages, or other communities. The gravity model could be applied with a 
structured-population (network) model, such as Kiskowskis model [43]. This model has been ap¬ 
plied to study community-based interventions [44]; it could provide a granular simulation of Ebola 
dynamics, with possible insights into community and differing levels of local intervention. 

During the course of the outbreak, parameters such as reporting rates, ETU availability, and 
intervention rates changed. Models that use a broader time-scale than our model see more benefits 
to incorporating time-variant parameter [22, 45]. In addition, transmission changes as behavior of 
infected and susceptible individuals changes [22]. Thus, models that incorporate behavioral changes 
could be successful in capturing that aspect of the epidemic. 

4.7 Conclusions 

We developed dynamic transmission models that account for spatial spread by considering the 
gravitational pull of larger populations and shorter distances in West Africa. The country-level 
model accurately captures epidemic dynamics and successfully forecasts cases and deaths for all 
three countries simultaneously. The district-level model captures the progression of EVD between 
and within districts in West Africa. Our models suggest that reduction of local transmission in 
Liberia and reduction of long-range transmission in Sierra Leone were the most effective interven¬ 
tions for the outbreak. The models also reveal differences in transmission levels between the three 
countries, as well as different at-risk population sizes. Our gravity spatial model framework for 
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EVD provides insight into the geographic spread of EVD in West Africa and evaluate the relative 
effectiveness of interventions on a large, heterogeneous spatial scale. Ultimately, gravity spatial 
models can be applied by public health officials to understand spatial spread of infectious diseases 
and guide interventions during disease epidemics. 
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5 Supplementary Information 

5.1 Model Details 

Country-level model. As discussed above, the model consists of S, E, R, I 2l F, and R compartments 
in each country. The model equations are given by: 

S n = —(Ci + £2 + £3) £>71 
En = (Ci + C2 + £3)871 — cnE n 
1 1 j> — O.E 7 1 q'n^ln Lgn-^lra 

ft-2 n — 'ynRn dl 2 n f2,n^2n 

En — SI 2 n S 2 F n 
Rn = C1 :ri 1 1 n T2,n^2n 
ICm = aE n 

E)C n = 6I 2n 

Cl — PljU^ln T ft2,nft2n T ft F.vEn 

C 2 — d n ,m(ft\.nftlrn T ft2.rift2rn T ftF,nErn) 

C3 = On,l(ftl,nhl + ft2,nhl + ftF, n Fl ) 

where 0 nm = K n ,ft nprr \ L . . 9 n 1 = n n ft np \ L , and n indicates each patch (Guinea, Sierra Leone, Liberia), 

\ a n,m) ’ \ a n,l) v ' 

and m and l the other two patches (with m l n). The risk of new infections from within each 
patch was represented by the £1 term for the home patch. TThe risk terms for EVD transmission 
from outside the patch, £2 and £ 3 , include 6, the gravity term. The distance exponent of the gravity 
term, t, is fixed at two. A range of values for 1 were tested, however we found that changes in t were 
compensated by changes in the fitted value of n. Thus, l was fixed for clarity, based on values in 
previous gravity models [28] . There are separate k for each country, so the rate of transmission into 
each country is different. This reflects the differing border porosities and rates of travel between 
countries. 


District-level model. Similarly, the district-level model equations for patch n are given by: 

S n = —(£1 + £,2)S n 
En = (£1 + £2 )S n - aE n 

Rn — OiE n 'fed In \ .<: l \ n 
l2n = 'Ycdln dl 2 n V2.c 12n 

En = dl 2 n S 2 F n 

Rn — C l,eft 1 n F 2 .C 12n 

ftC\n — OtE n 

E)C n = SI 2n 

£l — ftl,nft\n T ft2,nft2n T ftF,nE n 
63 

£2 = ^ ' ftl,nftlm T ft2,nft-2m T ftF,nE m 

m= 1 

ft - K P n P m 
(Sm)‘ 


( 2 ) 
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where Q n)U := 0 and c denotes the country to which patch n belongs. 

This model structure is similar to that of the country-level model. The main difference is that 
the long-range transmission is grouped into one £2 term, which sums the long-range transmission 
of every other patch into patch n. 

Model Parameters. 

The model parameters were estimated by simulating total cumulative cases, defined as the 
integral of the incidence ( aE n ) multiplied by k norm , a correction factor for underreporting and the 
fraction of the population at risk, among other factors [30]. 

Initial Conditions. 

Initial conditions for both models were determined based on the initial values of the data as 
follows: the number of susceptible persons in each patch was estimated based on the population 
size scaled by k norm . The total population of each patch was determined using data from the World 
Bank and national censuses [46-49]. The number of exposed persons was determined as twice the 
initial number of infected (ii and I 2 ) persons, based on the IZq for EVD, which has been estimated 
to be approximately 2 [16, 50]. The number of infected persons in II was determined based on the 
number of new cases in the previous nine days, based on the incubation period for EVD, and the 
number of infected persons in 12 was determined by subtracting the number of deaths within the 
next four days from the number of currently infected persons at the starting date [ ]. The number 
of infected persons from local versus long-range transmission was estimated based on the origin of 
the outbreak and the location of cases up until the starting point for the data: since the outbreak 
began in Guinea, the initial cases in Guinea were considered from local transmission, while the 
initial cases in Sierra Leone and Liberia were considered to be from long-range transmission from 
Guinea. The number of recovered persons was estimated based on the number of infected persons 
who did not die within the nine-day period that preceded the starting date. The number of persons 
in the F class was based on a fraction of the number of persons who died within the period before 
the starting date, estimated using the burial rate in Table 1. The initial conditions were all scaled 
to fractions of the population using the parameter k norm . 

We tested the country-level models ability to fit the outbreak when data was incorporated using 
a start date of March 30. In the simulations starting on March 30, local transmission in Liberia was 
turned off until May (when the first case appeared in Liberia), as the ODE framework used here 
cannot capture the stochastic process which leads to emergence of an outbreak in a new locale. The 
resulting fits and forecasts from these simulations were similar to those incorporating data from 
May onward (Supplemental Figures SI, S2). 

For the district-level model, parameters were determined using sampling from within plausible 
ranges and from the best-fits of the three-patch models. The 1 transmission rates were determined 
at the country level, reflecting the differing risks of transmission in the patches in each country. Each 
knorm was calculated to match final outbreak size, then adjusted by hand. The initial conditions 
were determined based on data from the WHO reports and from data compiled by the UN. The 
initial conditions were determined similarly to the country-level model, based on data from the 
WHO reports and from data compiled by the UN. 

The district capitals were considered to be the population centers for all districts except Montser- 
rado, Liberia, and Bonthe, Sierra Leone. Montserrado County is the location of the national capital 
and largest city, Monrovia, which was considered the population center. The capital of Bonthe is 
on a coastal island, so Yagoi, the closest city to Bonthe on the coast, was considered the population 
center. 

In the country-level model, interventions were simulated by reducing transmission parameters 
for either local or long-range transmission. Local transmission reductions were simulated by reduc¬ 
ing 1, the rate of local infections, in increments of one percent of the original value. Long-range 
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transmission reductions were simulated by reducing 2 and 3, the rates of infections from outside 
the patch, in one-percent increments. 

In the district-level model, interventions were simulated by eliminating the outbreak in one 
district (local transmission, long-range transmission into the district, and long-range transmission 
out of the district). The effect on the 62 surrounding patches was measured (difference in number of 
cases in those 62 patches in the normal model simulation versus number of cases in the 62 patches 
in the intervention simulation). In addition, the IAR score was determined as follows: 


IARScore = In 


Cases reduced in other districts 


- (3) 

70 of total outbreak occurring in intervention district 
Note that for districts with less than 0.05% of the outbreaks cases, the IAR score was defined 


to be 0. 
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5.2 Supplementary Figures and Tables 


Guinea Cases Sierra Leone Cases 




Liberia Cases 



Guinea Deaths Sierra Leone Deaths Liberia Deaths 





Figure SI: Cumulative cases and deaths for each country. Data depicted with black circles (from March 30, 
2014 to September 30, 2014) was used for model fitting to the data. Data depicted with red crosses (from 
October 1 to October 29) is displayed for validation of the model fits. The model projections are from March 
30 to October 31, 2014. The black circles represent the data used for fitting; the red crosses represent the 
data used for validation. The line represents the models overall best fit; the shaded regions represent the 
95% confidence interval. (Compare to main text Figure 2. 
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Guinea Sierra Leone 





Figure S2: Cases from local transmission and long-range transmission. Local transmission, from transmission 
within the country, is depicted in the color of the country in Figure 1. Long-range transmission, or cases 
from outside the country, is depicted in gray. Model projections are from March 30 to October 31. (Compare 
to main text Figure 3. 


20 












Guinea 






Liberia 


I” 

£ 

■ 

Gbarpolu, Liberia 

// 1 


/ s" 

//y u 

S” 

«20 


s= 1 

JE 

ze# - 



Isa • 

I* 

Jr ■: 




Figure S3: Cumulative cases and deaths for each district. Red lines represent cumulative cases (model) 
black lines represent cumulative deaths (model). Red and black circles represent data for cases and deaths 
respectively. Model projections are from March 30, 2014 to January 31, 2015. 











































































































































Figure S4: Comparison of all data points for cumulative cases and deaths in each of the 63 patches: com¬ 
parison of model values (x axis) and data (y axis). R 2 value of 0.8321. 
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Table SI. Parameter values and estimates for the country-level model. 


Parameter 

Definition 

Units 

Best Fit 

Range 

Sources 

a 

Transition rate from E to Ii 

days -1 

0.1059 

0.1-0.125 

[ . ] 

Pi G 

Transmission rate for I\ in Guinea 

persons -1 

days -1 

0.0950 

Initial: 0-0.5, Fitted: 0.0000 
to 0.0676 

fitted to data [ ] 

Pi SL 

Transmission rate for I\ in Sierra Leone 

persons -1 

days -1 

0.0504 

Initial: 0-0.5, Fitted: 0.0005 
to 0.1858 

fitted to data [ ] 

PlL 

Transmission rate for 1 1 in Liberia 

persons -1 

days -1 

0.1555 

Initial: 0-0.5, Fitted: 0.0727 
to 0.1656 

fitted to data [ ] 

PR 

Ratio of /1 transmission to I 2 and F 
transmission 

none 

2.4565 

1.5-3.0 

[5, 15] 

P2G,PfG 

Transmission rate for I 2 and F in 
Guinea 

persons -1 

days -1 

0.2335 

Calculated from Pi and pu 

See Pi, Pu 

P2SL , PfSL 

Transmission rate for I 2 and F in Sierra 
Leone 

persons -1 

days -1 

0.1238 

Calculated from Pi and Pu 

See Pi, P R 

P 2 L,Pfl 

Transmission rate for I 2 and F in Liberia 

persons -1 

days -1 

0.3820 

Calculated from Pi and 

See Pi, Pu 

Ai G 

Mortality for infected persons in Guinea 

none 

0.6643 

Initial: 0.5-0.9, Fitted: 0.3601 
to 1.4023 

fitted to data [ ] 

Al SL 

Mortality for infected persons in Sierra 
Leone 

none 

0.3910 

Initial: 0.5-0.9, Fitted: 0.2547 
to 0.6159 

fitted to data [ ] 

A L 

Mortality for infected persons in Liberia 

none 

0.6710 

Initial: 0.5-0.9, Fitted: 0.5453 
to 0.8236 

fitted to data [ ] 

^2 

Fraction of persons in I 2 that dies 

none 

0.9732 

0 .9-1.0 

[5, 7] 

5 

Death rate for persons in I 2 

days -1 

0.7597 

Calculated from A 2 and ko 2 

see A 2 and fco 2 

82 

Burial rate 

days - 

0.7425 

0.3333-1.0 

[29] 

ko 2 

1 /duration of I 2 

days -1 

0.7806 

0.3-0 .8 

[5, 32] 

k 2 i 

1 /duration of I\ 

days - 

0.1822 

0.1429 -0.2 

[ 1 ] 

i 

Distance exponent 

unitless 

2 


fixed 

Ki 

Between-patch coupling strength for 
Guinea 

km‘ 

persons -2 

10 -1017 

Initial: 10 - 8 -10 - 12 ,Fitted: 

10 -6 ' 21 - io -12 - 54 

fitted to data [ ] 

K 2 

Between-patch coupling strength for 
Sierra Leone 

km‘ 

persons -2 

10 -8 ' 27 

Initial: 10 - 8 -10 -12 , Fitted: 

IQ— 6.13_^q—21.80 

fitted to data [ ] 

K3 

Between-patch coupling strength for 
Liberia 

km‘ 

persons -2 

Iq-8.62 

lnitial:10 - 8 -10 -12 , Fitted: 

10 - 7 ' 89 - 10 -121 ° 

fitted to data [ ] 

tig 

Recovery rate of persons in I\ in Guinea 

days -1 

0.0578 

Calculated from Ai, A 2 , fei 

fitted to data [ ] 
(See Ai, A 2 , k 2 1 ) 

visl 

Recovery rate of persons in Ii in Sierra 
Leone 

days -1 

0.1090 

Calculated from Ai, A 2 , & 2 i 

fitted to data [ ] 
(See Ai, A 2 , fc 2 i) 
fitted to data 

til 

Recovery rate of persons in Ii in Liberia 

days - 

0.0566 

Calculated from Ai, A 2 , fei 

[l](See Ai, A 2 , 
k 2 1 ) 

t *2 

Recovery rate of persons in I 2 

days - 

0.0209 

Calculated from A 2 , ko 2 

[5, 7] 

7 g 

Rate and proportion of movement of 
persons from /1 to I 2 in Guinea 

days -1 

0.1243 

Calculated from Ai, A 2 , fei 

fitted to data [ ] 
(See Ai, A 2 , k 2 1 ) 

JSL 

Rate and proportion of movement of 
persons from /1 to I 2 in Sierra Leone 

days -1 

0.0732 

Calculated from Ai, A 2 , fei 

fitted to data 
[ ](See Ai, A 2 , 
k 2 i) 

7 L 

Rate and proportion of movement of 
persons from I\ to I\ in Liberia 

days -1 

0.1256 

Calculated from Ai, A 2 , fei 

fitted to data 
[l](See Ai, A 2 , 
k 2 i) 

See pi, p 2 ,pF 

Cl 

Transmission rate from S to Ii for lo¬ 
cally transmitted cases of EVD 

days -1 

Varies 

Calculated from Pi, P 2 ,Pf, 
h, h, F 

£ 2,6 

Transmission rate from S to Ji for long- 
range transmitted cases of EVD 

days -1 

Varies 

Calculated from Pi, p 2 ,pF, 
h, h, F 

See Pi, p 2 ,pF 

e 

Influence of long range cases on home 
patch 

unitless 

- 

Calculated from p, d, 1 , k 

See p, d, t, k 

Pu 

Population in patch n 

persons 

- 

Fixed for each country 

[ 3] 

dn,m 

distance between patch n and m 

km 

- 

fixed for each combination of 
patches 

[35] 

fcnorm 

Correction factor for reporting rate, 
population at risk, and other factors 

unitless 

0.0030 

0 . 001 - 0.1 

Sampled 
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Table S2. Parameter values for the district-level model. 


Parameter 

Definition 

Units 

Value 

Sources 


a 

Transition rate from E to I\ 

days” 1 

0.1059 

[5, 50] 



Pig 

Transmission rate for I\ in Guinea 

persons” 1 days” 1 

0.0950 

fitted 

[ ] 

to 

data 

Pi sl 

Transmission rate for I\ in Sierra 
Leone 

persons” 1 days” 1 

0.0504 

fitted 

[1] 

to 

data 

Pil 

Transmission rate for I\ in Liberia 

persons” 1 days” 1 

0.1555 

fitted 

[ ] 

to 

data 

Pr 

Ratio of I\ transmission to I2 and 
F transmission 

none 

2 

[5, 15] 



Ai G 

Mortality for infected persons in 
Guinea 

none 

0.6643 

fitted 

[1] 

to 

data 

&1SL 

Mortality for infected persons in 
Sierra Leone 

none 

0.3910 

fitted 

[1] 

to 

data 

A L 

Mortality for infected persons in 
Liberia 

none 

0.6710 

fitted 

[1] 

to 

data 

A2 

Fraction of persons in I2 that dies 

none 

0.97 

[5, 7] 



5 2 

Burial rate 

days” 1 

0.9 

[29] 



&02 

1/duration of I2 

days” 1 

0.8 

[5, 32] 



hi 

1/duration of I\ 

days” 1 

0.2 

[2, 32] 



7 G 

Rate and proportion of movement 
of persons from I\ to E in Guinea 

days” 1 

0.1370 

fitted to 
[1] (See 
A2, hi ) 

data 

Ai, 


Rate and proportion of movement 



fitted 

to 

data 

1SL 

of persons from I\ to I2 in Sierra 

days” 1 

0.0806 

[l](See Ai 

A2, 


Leone 



hi) 



7 L 

Rate and proportion of movement 
of persons from I\ to I\ in Liberia 

days” 1 

0.1384 

fitted to 
[l](See Ai 

& 2 l) 

data 

A2, 

i 

Distance exponent 

unitless 

2 

fixed 



hi 

Between-patch coupling strength 

km' - persons” 2 

10-S.5 

fitted 

[1] 

to 

data 

9 

Influence of long range cases on 
home patch 

unitless 

Calculated from 

p, d, i, k 

See p, 

d, t 

K 

Pn 

Population in patch n 

persons 

Fixed for each 
patch 

fixed for each 

[46] 



dn,m 

distance between patch n and m 

km 

combination of 
patches 

[35] 




Correction factor for reporting 






knorm 

rate, population at risk, and other 

unitless 

Varies by patch 

Fitted to data 


factors 
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Caption for Supplementary Animation: 

Video SI. A GIF representing the progression of Ebola according to the model from March 30, 
2014, to January 31, 2015. Color intensity of district represents number of cumulative cases in that 
district: darker color represents higher number of cases. 
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